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摘 要 : 以 塔里木 河 干流 为 研究 区 , 运 
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分 布 式 MIKE SHE 水 文 模 型 ,识别 了 日 值 尺 度 下 塔里木 河 干流 下 垫 面 各 参 


数 ,选取 干流 内 6 个 气象 站 (阿拉 和 尔 、 拜 城 . 库 车 ` 轮 台 库尔勒 EXT HL 0E) 1998—2007 年 实测 气象 数据 ,依据 干流 3 
个 水 文 站 (新 其 满 , 英 巴 扎 、 恰 拉 )1998 一 2007 年 实测 数据 对 模型 进行 识别 和 验证 ,并 结合 HadGEM 2 - ES 气候 模式 


中 气温 和 降水 数据 , 预 估 了 未 来 (2021 一 2050 4E) 气候 情景 下 塔里木 河 干流 径流 演变 特征 。 结 果 表明 :中 运用 


MIKE SHE 模型 模拟 塔里木 河 干流 日 径流 过 程 ,模型 效率 系数 均 达 0.72 以 上 ,表明 该 模型 适用 性 较 强 。@) 未 来 30 


a 塔里木 河 干流 多 年 
时 增加 趋势， 


] 平 均 气 温 均 较 基准 期 (1981 一 2004 年 ) 升 高 了 1.7 ~2.2 C , 且 夏 秋季 节 升 温 较 高 ;降水 整体 
春季 增幅 较 大 ,秋季 增幅 较 小 。@@) 未 来 气候 情景 下 ,各 水 文 站 年 均 径流 量 呈 减少 趋势 ， 


以 夏秋 季 


节 表 现 最 为 显著 (P<0.01)。 其 中 ,下 游 恰 拉 水 文 站 平均 径流 量 减少 比例 最 大 (5.04% ) ,上 游 新 其 满 水 文 站 减少 
比例 最 小 (0.6% ) ;各 气候 排放 情景 中 ,RCP2.6( 低 排放 情景 ) 下 各 水 文 站 减少 量 最 小 , RCP8.5 高 排放 情景 下 减少 
最 大 。 未 来 气候 情景 下 塔里木 河 干流 年 径流 呈 减 少 趋势 ,使 区 域 水 资源 调配 面临 更 大 困难 ,生态 安全 面临 更 大 威 


胁 , 同 时 也 将 加 大 区 域 人 地 关系 矛盾 。 
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径流 是 水 循环 的 基本 环节 ,同时 也 是 自然 地 理 
环境 中 最 为 活跃 的 要 素 。 日 径流 的 变化 模拟 与 分 析 
在 农业 灌溉 "工业 需 水 .水 资源 配置 与 评价 ”等 
方面 起 着 重要 作用 ,可 为 决策 者 提供 相关 参考 依据 。 
目前 ,国内 外 学 者 在 干旱 区 径流 模拟 方面 开展 了 一 
系列 的 研究 。 在 国外 , Abushandi 等 中 利用 HEC- 
HMS fil IHACRES 模型 分 别 对 Wadi Dhulied 干旱 区 
流域 径流 进行 了 模拟 ,表明 HEC-HMS 模拟 效果 优 
于 IHACRES;Worku 等 .中 利用 GIS 和 SWAT 模型 模 
JU f Beressa 流域 径流 和 沉积 物 对 土地 利用 的 啊 应 , 
结果 表明 土地 利用 方式 变化 显著 影响 径流 量 和 沉积 
物 量 。 黄 粤 等 "通过 运用 MIKE SHE 水 文 模型 模拟 
了 开 都 河流 域 的 日 径流 过 程 ,表明 MIKE SHE 模型 
能 够 应 用 于 干旱 区 ; 张 永 勇 等 ”基于 SWAT 模型 研 
究 了 石 羊 河中 上 游 流域 1988—2005 年 以 及 未 来 气 
候 情 景 下 出 山口 径流 量 的 变化 特征 ,探讨 了 未 来 气 
候 变 化 对 石 羊 河中 上 游 地 区 水 文 情 势 以 及 水 资源 的 
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影响 。 从 以 上 文献 及 水 文 模型 发 展 过 程 来 看 ,干旱 
区 径流 模拟 研究 主要 采用 的 模型 包括 SWAT 模型 、 
MIKE SHE 模型 TOPMODEL 模型 VIC 模型 等 。 本 
文选 取 塔 里 木 河 干流 为 研究 区 域 ,采用 MIKE SHE 
水 文 模型 研究 径流 变化 情况 。MIKE SHE 水 文 模型 
较 SWAT 等 其 他 模型 而 言 ,是 一 种 具有 物理 基础 的 
全 分 布 式 水 文 模型 ,能 够 用 来 模拟 全 部 的 水 文 过 程 ， 
同时 包含 灌溉 水 库 .水 工 建筑 物 等 优势 功能 。 塔 里 
木 河 干流 受 人 类 活动 影响 强烈 ,选取 MIKE SHE 模 
型 具有 明显 优势 ,塔里木 河和 干流 有 大 面积 的 灌区 , 
SWAT 模型 在 平原 灌区 应 用 时 存在 河 网 提取 、 人 工 
干扰 条 件 下 的 子 流域 和 水 文 单元 划分 以 及 作物 实际 
蒸 散发 量 计 算 等 问题 ?] ,因此 ,本 文采 用 MIKE SHE 
水 文 模型 进行 相关 研究 。 

气候 变化 作为 流域 水 资源 演 蔡 的 主导 ,在 区 域 
水 文 水 资源 调配 、 综 合 利 用 等 方面 具有 重要 影响 。 
伴随 着 全 球 变 暖 ,极端 气候 灾害 频 发 等 一 系列 气候 
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2 T 旱 区 


变化 现象 ,致使 气候 变化 对 水 文 水 资源 影响 范围 不 
靳 扩 增 ,影响 强度 日 益 加 大 ,国内 外 已 有 相关 学 者 对 
此 开展 了 相应 研究 。 在 国外 ,Christensen 4&7? 根据 
降 尺 度 模 拟 的 气候 情景 ,评估 了 气候 变化 对 克 罗 拉 
多 河流 域 水 文 水 资源 的 潜在 影响 ;Setegn 等 "利用 
全 球 气 候 模 式 输出 的 结果 研究 了 塔 纳 湖 盆 地 水 文 水 
资源 对 气候 变化 的 敏感 性 。 在 国内 , 张 建 云 等 ”1 指 
出 全 球 气候 变 暖 和 人 类 活动 直接 影响 了 水 循环 要 素 
的 时 空 分 布 特征 ,增加 了 极端 水 文 事件 发 生 的 概率 ; 
刘 昌 明 等 " 沁 曾 述 了 气候 变化 对 水 文 水 资源 的 主要 
影响 ,揭示 了 水 资源 系统 对 气候 变化 响应 的 主要 特 
征 及 研究 方法 ; 陈 亚 宁 等 "指出 气候 变 暖 引起 的 水 
资源 变化 ,将 会 使 西北 干旱 区 在 资源 开发 利用 过 程 
中 生态 维护 和 经 济 发 展 的 矛盾 更 加 突出 。 塔 里 木 河 
深 居 亚 欧 大 陆 内 部 , 属 典型 干旱 内 陆 河 ,而 塔里木 河 
干流 是 塔里木 河水 源 的 重要 耗 散 区 ,由 于 塔里木 河 
干流 水 资源 相对 紧缺 ,使 得 干流 沿岸 生态 十 分 脆 
弱 "”。 因 此 ,研究 未 来 气候 变化 下 塔里木 河 干流 径 
流 演变 特征 对 于 沿岸 社会 经 济 的 发 展 和 生态 环境 的 
保护 具有 重要 作用 。 本 文通 过 运用 MIKE SHE 分 布 
式 水 文 模型 模拟 塔里木 河 干流 日 径流 变化 过 程 , 探 
讨 气温 及 降水 将 如 何 影响 径流 过 程 以 及 各 河 段 水 量 
分 配 ,进而 预 估 了 未 来 30 a 塔里木 河 干流 各 河 段 水 
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量 调配 状况 , 旨 在 为 塔里木 河流 域 水 资源 配置 及 如 
何 应 对 气候 变化 对 径流 影响 提供 科学 依据 。 


1 研究 区 概况 


塔里木 河 干流 地 理 坐 标 为 81°13’ ~ 88?34' E, 
39°26' -41?32'N ,海拔 779 ~1 031 m ,面积 约 4.1x 
10* km  。 塔 里 木 河 干流 属 典 型 干旱 区 内 陆 河 ,自身 
不 产 流 , 历 史上 发 源 于 塔里木 河 盆地 周边 的 9 大 水 
系 都 曾 与 塔里木 河和 干流 有 着 水 力 联系 ,但 由 于 人 类 
活动 与 气候 变化 的 影响 , 现 仅 有 阿克苏 河和 田 河 、 
叶 尔 羌 河和 开 都 - 孔 省 河 与 塔里木 河 干流 有 水 力 联 
系 ,形成 了 目前 “四 源 一 干 ” 的 格局 “~ 干流 从 
肖 夹 克 到 台 特 玛 湖 全 长 1 321 km ,在 阿拉 尔 、 新 其 
满 、 英 巴 扎 、 乌 斯 满 恰 拉 分 别 建 有 水 文 观测 站 ,其 中 
阿拉 尔 水 文 站 多 年 平均 来 水 量 为 45. 82 x 10* m ,新 
其 满 水 文 站 多 年 平均 年 径流 量 为 37.62 x 10* m , 3 
巴 扎 水 文 站 为 28.76 x 10? m? , 恰 拉 水 文 站 为 6.78 x 
10* mm" 研究 区 气候 为 大 陆 性 暧 温带 极端 干旱 
沙漠 气候 ,年 平均 气温 在 10.6 ~11.5 C ,年 降水 量 
小 于 50 mm ,年 蒸发 量 2 300 ~3 000 mm ,光照 丰富 ， 
年 日 照 时 数 为 2 800 ~3 100 h ,无 霜 期 为 185 ~ 210 
d, 多 大 风 及 沙尘暴 天 气 , 区 内 植被 稀少 , 履 盖 度 
低 5 -2 o 
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图 1 塔里木 河 干流 边界 、 河 网 及 观测 站 点 分 布 


Fig.1 Boundary, river network and rainfall stations in main stream of Tarim river 


2 资料 与 方法 


2.1 气候 模式 情景 

第 5 次 耦合 模式 比较 计划 (Coupled Model Inter- 
comparison Project Phase 5, CMIPS ) 于 2013 年 发 布 
新 一 代 温 室 气 体 排放 情景 一 典型 浓度 路 径 ( Repre- 
sentative Concentration Pathways, RCP) ,RCPs 各 情景 


在 2100 年 选择 不 同 的 辐射 强迫 目标 ,主要 包括 
RCP 2.6( 低 排放 情景 )、RCP 4. 5/6. 0( 中 排放 情 
景 ) 及 RCP 8.5( 高 排放 情景 ) 共 4 种 排放 情景 。 本 
文采 用 英国 气象 局 哈 德 莱 中 心 HadGEM2 - ES 模式 
数据 ,假设 该 模式 对 当前 气候 模拟 的 优 劣 可 能 会 影 
响 对 未 来 气候 变化 的 预 估 结果 ,因此 要 检验 Had- 
GEM2 - ES 基准 期 的 气温 及 降水 数据 在 该 流域 降 尺 


4 期 


度 的 结果 " ,经 检验 得 出 气温 误差 较 小 ,降水 相对 
误差 较 大 。 为 了 提高 模式 数据 精度 ,本 文采 用 Delta 
方法 ”对 降水 数据 进行 校正 。 校 正 结果 如 表 1 所 
示 , 气 温 及 降水 误差 均 在 可 控 范围 ,表明 HadGEM2 
- ES 模式 数据 能 够 用 来 研究 气候 变化 对 塔里木 河 
干流 未 来 径流 的 影响 。 


表 1 HadGEM2 -ES 模式 数据 校正 结果 
Tab.1 Result of correction of HadGEM - ES data 


气象 站 点 气温 误差 /SC 降水 相对 误差 /% 
阿拉 尔 -0.01 1.8 
拜 城 0.02 1.2 
库 车 -0.07 0.4 
轮 台 -0.03 2.4 
库尔勒 0.08 3.1 
铁 干 里 克 0.02 0.8 


2.2 MIKE SHE 模型 

MIKE SHE 模型 是 全 分 布 式 且 具 有 物理 基础 的 
模型 系统 ,主要 用 于 模拟 坡 面 漫 流 .河流 和 湖泊 、 饱 
和 流 和 非 饱和 流 , 蒸 散发 .地 表 水 和 地 下 水 运动 及 其 
交换 等 方面 。MIKE SHE 模型 基于 单元 格 划 
分 流域 :在 水 平方 向 上 使 用 规则 格 网 离散 空间 ,垂直 
方向 上 利用 多 水 平 层 诠释 各 格 网 内 部 土壤 剖面 垂直 
差异 性 ,基于 数值 分 析 方 式 建 立 相 邻 格 网 要 素 间 以 
时 空 关系 ”1。 该 模型 主要 通过 6 个 独立 且 相 互联 
系 的 子 模块 构成 :截留 /蒸发 (ET) , 坡 面 漫 流 (OL)、 
河流 和 湖泊 计算 模拟 (0C)  AETRURT TEE DRE CUZ) .饱和 
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和 荒漠 6 类 ;灌溉 制度 和 灌溉 定额 参考 李 杰 友 "” 针 
对 干旱 区 水 资源 优化 配置 研究 的 结果 ,将 灌溉 面积 
按 水 文 站 点 所 属 区 位 划分 为 :阿拉 尔 一 新 其 满 、 新 其 
满 一 英 巴 扎 、 英 巴 扎 一 恰 拉 及 恰 拉 以 下 4 个 典型 河 
段 。 叶 面积 指数 参照 赵 成 义 、 韩 明 等 ”针对 塔 里 
木 河水 域 调配 的 研究 结果 ;土壤 类 型 数据 来 源 于 寒 
区 旱 区 科学 数据 中 心 (http ;// www. westgis. ac. cn), 
将 研究 区 土壤 类 型 划分 为 : 砂 壤 黏土 (sand. clay 
loam) XAK £ (silt loam ) 、 壤 精 土 (loam clay) 、 砂 
+ (sand) FAK + (clay loam)5 种 类 型 。 

本 人 研究 采用 FAO 推荐 的 Penman - Monteith 公 
式 , 对 干流 1998—2007 年 逐日 参考 作物 蒸 散发 量 进 
fip, 


900 
0. 408A(R, - G) *rg 275. 一 ev) 
A +r(1 +0. 34u, ) 


ET, = (1) 


UB LET, 为 参照 作物 蒸 散发 量 (mm * d^); R, 为 
作物 表面 上 的 净 辐 射 (MJ m? d ) ;6G 为 土壤 热 
通 量 (MJ m? * d);7 为 2m 高 处 日 平均 气温 
(C) ;u, 为 2m FS ABT XU (m. * s 7) ;e, 为 饱和 水 
汽 压 (kPa) ;e, 为 实际 水 汽 压 (kPa) ;A 为 饱和 水 汽 
压 曲 线 斜 率 ;r 为 湿度 计 和 常数 (kPa * 0C 7), 


表 2 塔里木 河 干流 MIKE SHE 模型 主要 输入 数据 
Tab.2 Data list of the MIKE SHE model of main 


stream of Tarim river 


带 流 (SZ) 、 融 雪 (SM) 模 块 , 且 每 个 子 模 块 均 对 应 一 
个 水 文 过 程 ,用 户 可 根据 研究 目的 自由 组 合并 搭建 
其 模块 。 
2.3 ”模型 数据 输入 

MIKE SHE 模型 需要 输入 数据 主要 包括 :DEM、 
逐日 气温 \ 降 水 数据 、 土 地 利用 数据 和 土壤 数据 等 。 
本 文 DEM 采用 比例 尺 为 1:250 000 的 栅 格 数据 (ht- 
tp://srtm. csi. cgiar. org/SELECTION/inputCoord. 
asp) 。 气 温和 降水 基准 数据 来 源 于 国家 气象 信息 中 
心 (http :// data. cma. en/ ) ,其 中 ,以 干流 及 周边 的 6 
个 气象 站 点 构建 泰 森 多 边 形 ,用 以 代表 整个 干流 降 
水 ,并 采用 联合 国 粮农 组 织 推 荐 的 Penman - 
Monteith 公式 计算 干流 参考 蒸 散发 值 。 土 地 利用 数 
据 来 源 于 中 国 科 学 院 资源 环境 科学 数据 中 心 ht- 
tp://www. resdc. cn) ,其 中 ,将 2005 年 划分 的 21 类 
土地 利用 数据 合并 为 农田 、 林 地 、 草 地 水 体 、 居 民 点 


数据 类 型 数据 来 源 属性 
地 形 1: 250 000 DEM 1 km x 1 km 格 网 
河 网 水 系 中 国 地 区 三 级 以 上 河流 。 * .shp 文件 
土地 利 nn 1 km x 1 km 格 网 
土壤 数据 FAO 土壤 图 1 km x 1 km 格 网 
降水 /气温 数据 ”国家 气象 信息 中 心 6 站 
径流 数据 塔里木 河流 域 管理 局 4 站 


2.4 模型 参数 率 定 
采用 1998—2007 年 塔里木 河 干流 日 径流 数据 
进行 模拟 ,模型 率 定期 为 1998 一 2003 年 ,模型 验证 
期 为 2004 一 2007 年 。 采 用 相对 误差 (RE) 、 相 关系 
数 (CORR) 和 纳什 系数 (NS ) 等 指标 评价 模型 适用 
pps ` 
S-0 


RE = 5 x 100% 


(2) 
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2 塔里木 河 干流 土壤 类 型 和 土地 利用 


Fig.2 Soil types map and landuse map of main stream of Tarim river 


Y (0, -万 (5 - 8) 
CORR - E - (3) 
> (0; -0) Y (S; - 5)? 


i=l 


0; - S) 
mand g 
2 (enr 


式 中 :0, 表示 第 CX Fill m «s S0 表示 实 
测 径流 量 平 均值 ,m «87:5, 表示 第 i 天 模拟 径流 
量 ,m :si 表示 模拟 径流 量 平均 值 ,mm s^, 


3 结果 


3.1 模型 识别 与 验证 

参考 Refsgaard!” ?" 提出 的 MIKE SHE 模型 参 
数 化 原则 ,将 参数 划分 为 两 类 :中 ) 实测 数据 或 台 感 
数据 推算 所 得 ,该 类 数据 是 基于 对 现实 情况 的 真实 
描述 ,存在 一 定 的 不 确定 性 , 故 在 目前 数据 可 获取 程 
度 无 法 提高 情况 下 ,将 降水 ,参考 蔡 散 发 .气温 、 土 地 
利用 及 植被 等 参数 不 予 率 定 ; D 非 测量 和 计算 所 
得 ,包括 河道 的 Manning PER .河床 透水 系数 、 植 被 
参数 .土壤 与 含水 层 参数 等 “” 。 本 文 结 合 以 上 参数 
率 定 原则 ,将 塔里木 河 干流 各 水 文 站 1998 一 2003 年 


实测 日 径流 数据 用 于 率 定 模型 ,并 按照 从 干流 上 段 至 

下 段 依次 率 定 的 顺序 进行 判定 ,其 结果 如 表 3 所 示 。 
由 表 4 可 知 , 率 定 期 各 水 文 站 相关 系数 均 二 0. 

87 ,纳什 系数 均 >0.72 ,相对 误差 均 < 70.07 ;验证 


表 3 参数 率 定 结果 
Tab.3 Result of values of the calibrated parameters 


模块 参数 河 段 率 定 结果 
坡 面 流 曼 宁 系 数 干流 0 
河道 流 ”河道 糙 率 系数 干流 30 

河床 透水 系数 。 阿拉 尔 - 新 其 满 。 1.1x10-° 

新 其 满 - 英 巴 扎  7.0x1077 

英 巴 扎 =- 恰 拉 7.5x1077 

恰 拉 - 台 特 玛 湖 7.5x1077 

饱和 流 渗透 系数 干流 6.8 x10-5 

给 水 度 干流 0. 125 
表 4 模型 结果 


Tab.4 Results of MIKE SHE model 


Sue d " TAS 
阶段 水 文中 p a Dm 
率 定 期 新 其 满 0.97 0.93 -0.07 
英 巴 扎 0.91 0.73 -0.36 
恰 拉 0.87 0.72 -0.11 
验证 期 新 其 满 0.97 0.94 -0.02 
英 巴 扎 0.94 0. 84 -0.21 
Tg du 0. 89 0. 79 -0.18 
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图 3 


日 期 


塔里木 河 干 流 各 水 文 站 日 径流 模拟 


Fig.3 Simulation of daily runoff of every station of main stream of Tarim river 


Hio K CX ROS ORO] zm:0. 89 , ZEE ROS] m0. 79 , 
模拟 与 实测 的 径流 曲线 拟 合 性 较 好 ,表明 通过 率 定 
参数 ,模型 能 较为 准确 地 模拟 塔里木 河 干流 日 径流 
变化 趋势 。 
3.2 未 来 气候 情景 下 径流 变化 预 估 
本 研究 将 前 述 的 塔里木 河 干流 及 周边 6 个 气象 
站 点 经 校正 后 的 HadGEM2 - ES 模式 数据 计算 其 平 
iis CR MEN c es 
的 总 体 趋势 。 根 据 图 4 中 气温 相对 变化 可 知 ， 
2021—2050 年 塔里木 河 干流 在 不 同 排放 情景 下 多 
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3 上 
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2 
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1 
0.5 
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图 4 2021—2050 年 塔里木 河 干流 多 年 


年 月 平均 气温 较 基 准 期 均 呈 上 升 趋势 ,升温 幅度 差 
异 显 著 ,集中 于 1 ~2 %C 间 。3 种 排放 情景 下 ， 
RCP8.5 情景 升温 幅度 明显 高 于 其 他 2 种 。 从 各 月 
升温 情况 来 看 ,夏秋 季节 增 温 较 大 ， 升温 最 大 为 1 月 
份 ,升温 幅度 为 1.5 ~3 C. Z% E ,2021—2050 年 
各 排放 情景 较 1981 一 2004 年 基准 期 多 年 月 平均 气 
温 升 高 趋势 明显 ,夏秋 季节 增 温 较 大 , 春 冬季 节 增 温 
较 小 ,整体 增 温 幅度 在 1.7 ~2.2% 间 。 由 图 4 可 以 
看 出 ,塔里木 河 干流 不 同情 景 下 2021—2050 年 多 年 
月 平均 降水 较 基 准 期 1981 一 2004 年 多 年 月 平均 降 


50 pf EIRCP2.6 ORCP4.5 
40 F 


RCP8.5 


降水 量 相 对 变化 % 


1:28: A e 07 8 9 TO ll 12 
月 份 


月 平均 气温 降水 较 基准 期 变化 量 


Fig.4 Change of multi — month mean temperature and precipitation from 2021to 2050 relative to the base period 


水 呈现 春季 增幅 较 大 ,秋季 增加 较 小 趋势 。 从 各 月 
来 看 ,2、5、8、9、11 月 有 部 分 气候 情景 的 降水 量 减 
少 ,其 他 各 月 各 气候 情景 降水 量 都 是 增加 的 。 

由 于 塔里木 河 干流 是 纯 耗 散 性 河流 ,自身 不 产 
流 ,在 气候 模式 下 需要 计算 出 阿克苏 河 、 和 田 河 、 叶 
尔 羌 河 总 径流 量 才能 作为 气候 模式 下 阿拉 尔 站 的 输 
Adi ,模式 情景 工作 量 较 大 且 较 难 控制 。 然 而 ,本 研 
究 区 阿拉 尔 站 1957 一 2000 年 的 多 年 平均 年 径流 量 
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图 5 塔里木 河 干流 2021—2050 年 各 水 文 站 多 年 
平均 径流 较 基 准 期 2003 年 相对 变化 


Fig.5 Change of multi — annual monthly mean runoff amount 


of every station from 2021 to 2050 relative to the base period 


35 4 


H 45. 82. x 10° mi ,阿拉 尔 站 2003 年 年 径流 量 达 44. 
98 x 10* m? ,这 与 多 年 平均 径流 量 较为 接近 ,因此 ， 
本 文选 取 阿拉 尔 站 2003 年 径流 量 为 基准 期 输入 模 
型 ,同时 将 CMIP 5 中 HadGEM2 - ES 情景 模拟 的 6 
个 气象 站 的 RCP2. 6, RCP4. 5 和 RCP8.5 下 降水 和 
最 高 最 低 平 均 气 温 数 据 作为 模型 的 输入 值 ,进而 驱 
动 MIKE SHE 模型 输出 塔里木 河 干流 各 水 文 站 径流 
量 。 由 表 5 和 图 5 可 知 , 未 来 气候 情景 下 ,新 其 满 、 
英 巴 扎 和 恰 拉 3 个 水 文 站 年 径流 量 较 基准 期 2003 
年 呈现 出 减少 趋势 ,其 中 恰 拉 水 文 站 减少 比例 最 大 ， 
达 5.04% ,新 其 满 水 文 站 减少 比例 最 小 , 达 0. 6% 。 
从 各 气候 排放 情景 可 知 ,在 RCP2.6 情景 下 各 水 文 
站 减少 量 最 小 ,在 RCP8. 5 情景 下 各 站 减少 量 最 大 。 
由 于 春 冬 季节 气温 增加 ,蒸发 增 大 ,同时 春 冬 季 径 流 
量 较 小 ,使 得 各 水 文 站 径流 量 春 冬季 节 减 少 比 例 较 
为 显著 。 夏 秋季 节 增 温 明 显 , 莱 发 加 强 , 各 水 文 站 径 
流量 各 情景 下 减少 量 相对 较 大 ,但 由 于 夏秋 季节 径 
流量 较 大 ,使 得 相对 变化 较 春 冬 季 市 变化 微弱 。 对 
于 降水 而 言 , 虽 表 现 为 增加 趋势 ,但 塔里木 河 干流 降 
水 基本 不 产生 地 表 径 流 , 是 水 资源 的 主要 耗 散 和 消 
耗 区 。 因 此 ,未 来 气候 情景 下 气温 对 各 水 文 站 径流 
量 相 对 变化 影响 较为 显著 。 


表 5 2021 -2050 年 塔里木 河 干流 各 水 文 站 多 年 径流 量 
较 基 准 期 2003 年 相对 变化 ( %) 
Tab.5 Change multi - annual mean runoff amount 


of every station from 2021 to 2050 relative 


to the base period / 96 
各 气候 情景 新 其 满 英 巴 扎 恰 拉 
RCP2.6 -0.60 -1.44 -4.54 
RCP4.5 -0.62 -1.54 -4.82 
RCP8.5 -0.65 -1.64 -5.04 


4 结论 与 讨论 


本 文 基 于 MIKE SHE 模型 探究 了 塔里木 河 干流 
日 径流 变化 特征 ,并 结合 HadGEM2 - ES 模式 数据 
对 径流 响应 进行 了 模拟 分 析 ,结果 表明 : 

(1) 通过 对 塔里木 河 干流 径流 模拟 结果 与 实测 
值 的 验证 , 率 定期 塔里木 河 干流 3 个 水 文 站 日 径流 
相关 系数 均 在 0. 87 以 上 ,纳什 系数 均 在 0.72 以 上 ; 
验证 期 各 水 文 站 日 径流 相关 系数 均 在 0.89 以 上 , 纳 
什 系数 均 在 0.79 以 上 ,表明 ,MIKE SHE 水 文 模型 
在 干流 模拟 的 可 行 性 较 好 。 这 与 田 开 迪 等 “利用 
MIKE SHE 模型 对 灞 河流 域 径 流 模拟 的 结果 相 一 
致 , 均 表明 MIKE SHE 模型 对 流域 模拟 具有 较 好 的 


ChinaXiv 合 作 期 刊 


4 期 班 春 广 等 :未 来 气候 情景 下 塔里木 河 干流 日 径流 过 程 模拟 研究 7 


适用 性 ,但 本 文 重点 针对 日 径流 模拟 且 模 拟 效果 较 
好 ,而 田 开 迪 等 * 的 研究 以 年 径流 为 基准 ,并 未 探 
究 月 径流 和 日 径流 模拟 效果 ,其 研究 区 也 以 半 湿 泣 
气候 为 主 。 因 此 ,针对 MIKE SHE 模型 对 干旱 内 陆 
河流 域 的 适用 性 ,还 有 待 于 下 一 步 从 月 径流 .年 径流 
等 多 尺度 综合 探究 。 

(2) 各 情景 模式 下 多 年 月 平均 气温 较 基准 期 呈 
著 增 加 趋势 ,上 且 以 夏秋 季节 增 温 较 大 ,3 种 情景 升 
温 幅 度 在 1.7 ~2.2 % ;降水 在 各 水 文 站 各 情景 下 呈 
整体 增加 趋势 , 旦 春季 增幅 较 大 ,秋季 增幅 较 小 , 表 
明 未 来 气候 可 能 将 趋 于 暧 一 湿 状 态 。 这 与 陶 辉 
等 通过 采用 多 模式 加 权 平 均 数据 对 塔里木 河流 
域 未 来 气温 变化 的 研究 结果 相 一 致 , 均 表明 未 来 气 
候 基 于 3 种 排放 情景 下 塔里木 河流 域 气温 均 呈 上 升 
趋势 ,升温 幅度 均 在 2 CC 左右。 与 此 同时 ,本 结果 还 
与 罗 敏 等 “针对 和 田 河流 域 未 来 不 同 排放 情景 
春秋 、 冬 三 季 降 水 均 增加 ,夏季 降水 减少 ,温度 升 高 
且 夏 季 增 温 最 大 的 研究 结果 一 致 。 年 际 尺度 下 , 塔 
里 木 河 气温 和 降水 均 呈 不 同 程度 的 增加 趋势 ,这 与 
Liu 457" fidi 18 种 全 球 气 候 模 式 得 到 的 多 模式 集 
合 平均 数据 ,对 喀 喇 昆仑 山地 区 气温 降水 的 研究 结 
果 相 一 致 , 均 表明 在 RCP4.5 和 RCP8.5 情景 下 气温 
和 降水 呈现 不 同 程度 的 增加 趋势 。 

(3) 3 个 水 文 站 径流 量 均 有 所 减少 ,其 中 恰 拉 
水 文 站 减少 比例 最 大 , 达 5. 04% ,新 其 满 水 文 站 减 
少 比 例 最 小 , 达 0.6% 。 由 于 气温 升 高 ,整个 干流 蒸 
发 加 强 ,使 得 各 水 文 站 径流 量 呈 现 出 不 同 程度 的 减 
少 趋势 。 各 气候 排放 情景 中 , RCP2. 6 情景 下 各 站 
减少 量 最 小 ,RCP8. 5 情景 下 则 呈现 为 最 大 。 未 来 
塔里木 河 干流 由 于 气温 升 高 ,蒸发 加 强 ,致使 各 水 文 
站 径流 量 将 趋 于 减少 态势 。 这 与 罗 敏 等 “针对 和 
田 河 流域 未 来 不 同 排放 情景 下 ,7 一 8 月 径流 量 呈 不 
断 减 少 的 结果 相 一 致 , 均 表 明和 干旱 内 陆 河 干流 区 未 
3€ 30 a 径流 量 呈 减少 趋势 。 这 一 结果 一 方面 将 加 
剧 该 区 域 水 资源 短缺 ,给 塔里木 河 干流 水 资源 优化 
配置 带 来 更 大 挑战 ; 男 一 方面 将 使 生态 安全 形势 更 
加 严峻 ,致使 生态 稳定 性 中 源 性 较 弱 的 塔里木 河流 
域 面临 更 大 的 危机 。 

(4) 本 文选 择 降 尺度 的 方法 探究 塔里木 河 干流 
径流 变化 过 程 , 然 而 降 尺 度 方法 的 选择 模式 数据 校 
正方 法 的 选择 .单一 模型 及 单一 气候 模式 为 依据 , 研 
究 气候 变化 对 水 文 水 资源 和 环境 的 影响 都 存在 较 大 
的 不 确定 性 5 ,致使 模拟 预测 结果 存在 一 定 的 误 
差 。 因 此 ,提高 区 域 气候 变化 情景 预测 精度 , 综 


zi 


合 多 模式 多 模型 研究 方法 ,探究 气候 变化 对 水 文 水 
资源 的 影响 ,可 大 大 降低 预测 的 不 确定 性 。(5) 本 
文 对 于 未 来 30 a 塔里木 河 干流 各 河 段 径流 演变 特 
征 的 研究 , 仅 将 塔里木 河 干流 多 年 平均 年 径流 量 输 
人 模型 ,未 能 考虑 丰 水 年 . 枯 水 年 各 河 段 径流 变化 情 
况 ,下 一 步 将 利用 构建 的 塔里木 河 干流 径流 模型 , 结 
合 以 上 数据 进行 丰 水 年 、 枯 水 年 各 径流 演变 方面 的 
探究 ,进一步 完善 塔里木 河 干流 水 资源 合理 配置 及 
水 资源 评价 。 

致谢 :感谢 审 稿 专家 对 本 文 修改 所 提出 的 建设 
性 意见 和 编辑 人 员 耐 心细 致 的 指导 。 
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Simulation Research of Daily Runoff Process in main stream of Tarim 
River under Scenario of Future Climate 
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Abstract: Taking main stream of Tarim River as research area, applying hydrological model of distributed MIKE 
SHE , various underlying surface parameters of main stream of Tarim River under daily value dimensions were identi- 
fied,six meteorological stations ( Alar, Baicheng , Kuche , Luntai, Kuerle, Tieganlike) observed datas from 1998 — 
2007 were used ,the observed discharge at three hydrometric stations ( Xin Qiman , Ying Bazha ,and Chara) observed 
datas from 1998 —2007 were used for the calibration and validation , and combining the air temperature and rainfall 
data in climate model of HadGEM 2 - ES ,so as to estimate the future (2021 to 2050) runoff evolvement features of 
main stream of Tarim River under climate scenarios. The results indicated that :1. The model of MIKE SHE was a- 
dopted to simulate runoff process of main stream of Tarim River , with model efficiency coefficients over 0. 72 , which 
showed that this model had a strong applicability. 2. In the future 30 years , the monthly average temperature of main 
stream of Tarim River will be increased by 1. 7 -2.2°C compared with that during the baseline period (1981 to 
2004) ,with obvious temperature rise in summer and autumn. The rainfall overall showed an increasing trend with a 
higher increase in spring and lower increase in autumn. 3. Under the future climate scenarios , annual runoff volume 
of various hydrometric stations shows decreasing tendency , with obvious decrease in summer and autumn (P < 0. 
01). The downstream Chara hydrometric station has the biggest decrease rate (5. 0496 ) , while the upstream Xin Qi- 
man hydrometric station presents the opposite (0.696 ). In emission scenarios under various climates ,as opposite to 
various hydrometric stations with RCP 8.5 ,that with RCP 2.6 has the smallest decrease rate. The annual runoff of 
main stream of Tarim River under the future climate scenarios showed decreasing tendency , causing greater difficul- 
ties in allocating the regional water resources and greater threat to the ecological security. In the meanwhile ,the con- 
tradictions of regional man - earth relationship will also be aggravated. 


Key words: Daily runoff model; Main stream of Tarim River; Climate change; MIKE SHE 


